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In this paper we study the isocurvature mode contribution to the cosmic microwave background 
anisotropics and the large scale structure power spectrum, for a two-field model of inflation proposed 
by Linde and Mukhanov. We provide constraints on the parameters of the model by comparing its 
predictions with observations of the microwave background anisotropics, large scale structure data 
on the galaxy power spectrum, and the number density of nearby galaxy clusters. We find that 
such models are consistent with observations for a narrow range of parameters. As our main result, 
we find that only a very small isocurvature component is allowed, a < 0.006, for any underlying 
Friedmann model. Furthermore, we give the expected resolution with which the model parameters 
will be determined from future satellite missions like MAP and Planck, for a fiducial flat ACDM 
model. We find that Planck mission will be able to detect such small contributions, especially if 
polarization information is included. The isocurvature spectral index m so will also be determined 
with better than 8% precision. 
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I. INTRODUCTION 

The future satellite experiments MAP Q| and Planck 
|^[ open an exciting era for cosmologists and for particle 
physicists. The high resolution and sensitivity of these 
experiments will allow such a precise determination of 
the cosmic microwave background (CMB) power spec- 
trum of temperature and polarization anisotropics that 
it will soon be conceivable to test different cosmologi- 
cal models with great accuracy ||. Meanwhile, galaxy 
surveys, like 2dF and Sloan Digital Sky Survey (SDSS), 
aimed at measuring several hundred of thousand redshifts 
will provide a map of the Universe with unprecedented 
precision and extension. Any viable cosmological model 
must produce reasonable fits to both the observed CMB 
and large scale structure (LSS) power spectra, and must 
be tested on the basis of all available data. 

Within the context of inflationary scenarios, Gaussian 
adiabatic fluctuations are often assumed as a standard 
prediction. However, besides the usual adiabatic fluc- 
tuations, other independent modes may be present, e.g. 
isocurvature fluctuations 0-0, and in many cases they 
are non-Gaussian Q. In particular, the isocurvature 
mode is an entropy perturbation, characterized by an 
appropriate balance of the fluctuations in the different 
components, such that the spatial curvature remains un- 
perturbed. There is an attractive model that has been 



proposed some time ago that considers pure isocurva- 
ture perturbations in the baryon component but un- 
fortunately seems to be ruled out by present observa- 
tions [Plpl ■ Most of the recent models are actually cold 
dark matter (CDM) isocurvature models fl^ |, but there 
is also an ingenious neutrino isocurvature model jll|. 

Far from being academic, the reason for considering 
also isocurvature fluctuations resides in the fact that 
many different inflationary models, with more than one 
scalar field, predict the formation of significant isocurva- 
ture fluctuations during the inflationary era [§,|l4| . As for 
their predictions on the density perturbation power spec- 
tra, they may differ among themselves in having different 
amplitudes and tilts of isocurvature and adiabatic spec- 
tra, as well as for the statistical nature of both modes. 
Many inflationary models predict that both the isocurva- 
ture and adiabatic fluctuations have a nearly scale invari- 
ant (Harrison-Zeldovich) spectrum with Gaussian statis- 
tics. These models have been tested against LSS in the 
recent literature (l^,[l(| , which showed that only a small 
fraction of isocurvature component seems to be allowed 
by present observations. 

In this paper we focus on a particular inflationary sce- 
nario proposed by Linde and Mukhanov fi~7| , in which 
none of the two conditions mentioned above is neces- 
sarily satisfied. In fact, here the isocurvature fluctua- 
tions are non-Gaussian, more specifically, they are x 2 
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distributed]^] and their spectrum is a power law, with 
spectral index ri; so > 1. We test this model with the 
present observations on CMB and LSS, and we also make 
predictions on how well the future satellite experiments 
will be able to measure the relevant model parameters. 

In previous papers |ll||l6|,|l8| the mixed spectra of 
adiabatic plus isocurvature modes were assumed to be 
independent, Gaussian and approximately scale invari- 
ant. Some 1 15 1(| used the large scale structura data, 
together with the COBE-DMR normalization, to con- 
strain the models, while others jl8| only used the CMB 
anisotropics, with or without polarization. In our pa- 
per, we have combined both CMB and LSS observations, 
for a non scale invariant spectrum of perturbations, and 
included also the gravitational wave contribution. 

In section |n| we review how isocurvature modes 
may arise from inflation, and we describe the Linde- 
Mukhanov model which inspired this work. In section III 
we introduce the power spectra for matter and radiation 
in the mixed (adiabatic + isocurvature) case, showing the 
effect of the amplitude and tilt of the different isocurva- 
ture contributions. In section |iy| we show the compar- 
ison of the mixed spectra with the available data: we 
constrain the parameter range considering a joint analy- 
sis of both CMB and LSS, and we estimate with which 
precision the future satellite experiments will determine 
the relevant parameters. Finally, section ^ is dedicated 
to a general discussion of the results. 



II. ISOCURVATURE MODES FROM INFLATION 

Isocurvature perturbations are generated during infla- 
tion whenever there is more than one scalar field present. 
They correspond to entropy perturbations that do not 
perturb the metric, and thus the spatial curvature. They 
typically arise when one of the fields is fixed by its poten- 
tial during inflation, the infiaton energy later decays into 
relativistic particles and redshifts away, while the other 
field's energy becomes the dominant contribution. De- 
pending on model parameters, the relative contribution 
of adiabatic to isocurvature perturbations may be notice- 
able in the microwave background anisotropies and large 
scale structure, and they may have in principle very dif- 
ferent spectral tilts, e.g. blue (n > 1) or red (n < 1). 
Furthermore, depending on the evolution during infla- 
tion, the statistics of the different components (adiabatic 
and isocurvature) could be very different, e.g. Gaussian 
and x 2 distributed, respectively. Such a complicated phe- 
nomenology requires a detail analysis in order to confirm 
whether a particular model is ruled out by observations. 



1 However, slight variants of the model can also give a Gaus- 
sian isocurvature mode j^j. 



In this paper we will concentrate in a particular re- 
alization of a mixed adiabatic and isocurvature model 
proposed recently by Linde and Mukhanov. 



A. The Linde— Mukhanov model 

The model of Ref. |0] is probably the simplest one can 
think of that produces isocurvature perturbations during 
inflation. It has two coupled massive scalar fields de- 
scribed by the scalar potential, 
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In principle, inflation could occur along either of the val- 
leys at <f) = or a = 0, depending on initial conditions. In 
a chaotic inflation approach one expects that the fields 
will start at very large values, <p, a 3> Mp, where the 
coupling term g 2 a 2 (j> 2 dominates. Let us suppose that 
initially one of the fields has a larger value, say \c/>\ > \a\ 
and thus the field a rapidly settles at a = 0. Then infla- 
tion occurs along the a = valley, with energy density 



^M 2 (j> 2 and a Hubble constant 



H z = 
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During inflation the mass of the a field becomes 



vH 2 
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where v — 3g 2 Mp/4:irM 2 is a constant. Typically, during 
inflation the second term dominates and thus the model 
gives a mass term of the a field proportional to the rate 
of expansion. 

The quasi-de-Sitter evolution during inflation provides 
a neat way to generate metric perturbations from quan- 
tum fluctuations. Those of the infiaton will give rise to 
adiabatic density perturbations, since the energy den- 
sity during inflation is proportional to the infiaton field 
fluctuations, Sp ~ V'((f>)S(j). On the other hand, quan- 
tum fluctuations of the a field will not generate curvature 
perturbations since the inflationary trajectory lies along 
(a) = 0. Nevertheless, after inflation the energy den- 
sity of the a field may come to dominate the evolution 
of the universe (e.g. as a cold dark matter component) 
and its fluctuations would then contribute as isocurva- 
ture perturbations 0]. Let us compute the amplitude of 
those perturbations. For a massive field with m 2 <C H 2 
during inflation, the amplitude of the long wavelength 
perturbation of the a field at the end of inflation is given 
by 
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and the average perturbations of energy density in the a 
field, 8p a = rh 2 (Sa) 2 /2 can be estimated as 
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which corresponds to a "blue" spectral index 
n iso w 1 + 4m 2 /3-ff 2 ~1 + 4^/3 > 1 . 



(5) 



(6) 



The small ratio fh 2 /Mp ensures that the a field does not 
contribute initially to the perturbations of the metric, 
i.e. it generates isocurvature perturbations, which much 
later could end up dominating, as mentioned above. The 
detailed evolution is very model dependent jl7) . We will 
assume, following Linde and Mukhanov, that the correc- 
tion vH 2 to the mass of the a field disappears soon after 
inflation, when the inflaton field <$> decays into relativistic 
particles while the a field remains stable or decays very 
late, and thus its energy density (in coherent oscillations 
of the field) may dominate in the form of cold dark mat- 
ter today. Under these assumptions one can estimate the 
corresponding density contrast J17| 



Sj ^-C(k) (±- 



2i//3 



(7) 



for k <C mm[My/m/M, M exp(— l/2v)], where 
C{k) 
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This isocurvature perturbation has non-Gaussian statis- 
tics, in fact a \ 2 distribution, because it arises from the 
square of a Gaussian field [see Eq. (j|)]. f\ 

On the other hand, the adiabatic density perturba- 
tions generated by the inflaton field <f> during inflation 
contribute with a Gaussian spectrum with amplitude 
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where Nk = A^or = 65 is the number of e-folds before 
the end of inflation when the mode with wavenumber k 
corresponding to our present horizon crossed the Hubble 
scale during inflation, with the spectral index 



"-ad 



1 



2 



(10) 



There is enough freedom in the model parameters 
(u, M, m) to have the Gaussian adiabatic perturbations 



One could avoid having non-Gaussian statistics for the 
isocurvature component by giving the a field a vacuum ex- 
pectation value, (a) 7^ 0. In this way, the main isocurva- 
ture contribution would come from the term 8p a — m 2 (a)5a, 
which is Gaussian distributed. 



dominate the spectrum on large scales, while the non- 
Gaussian isocurvature perturbations dominate on small 
scales. In practice, the latter will never dominate at 
scales of cosmological interest since, as we shall illustrate 
in Sections IV and V, the relative contribution of the 
isocurvature component is constrained by present obser- 
vations to be very small. 

We will thus consider a mixed model of adiabatic and 
isocurvature perturbations arising from inflation and con- 
tributing simultaneously to the temperature anisotropics 
of the microwave background and the power spectrum of 
large scale structure. The relative tilts and amplitudes 
of the two spectra will allow us to compare with present 
observations and determine the likelihood functions for a 
given set of model parameters. 

We will use the CMBFAST code HQ to normalize 
the model to COBE data [^l|-^3|, compute the C; and 
power spectra P{k), and then compare with observations. 

Therefore, the post-inflationary power spectra for the 
independent adiabatic and isocurvature components read 



(11) 
(12) 



P ad (k) = 25 x 10- iu ( — 

. "'hor 



Piso(fc) - dio x 10- iU ( — 

Kgal 



where dio =0(1 — 20) is an arbitrary normalization fac- 
tor, to be fixed by observations. We will discuss in the 
next Section the relative normalization of the two com- 
ponents for a mixed power spectrum. 



III. MATTER AND RADIATION POWER 
SPECTRA 

In order to compare the theory with the data, we com- 
pute the power spectra of matter and radiation at the 
present epoch. For this purpose, we solve Boltzmann 
equations using CMBFAST, starting from a very early 
epoch [z ~ 10 7 ) up to the present time. 

In the presence of only one mode of fluctuations the 
actual matter power spectrum is: 



P(k) = T 2 (k)P in (k), 



(13) 



where Pi n (k) denotes the initial power spectrum and the 
transfer function T(k) contains the information on how 
the spectrum is modified through the evolution. T(k) 
is typically of order 1 on large scales (small k) and it 
becomes less than 1 at smaller scales (larger k), the ac- 
tual shape depending on the specific mode considered 
(i.e., whether adiabatic or isocurvature), the cosmologi- 
cal parameters and on the nature of the dark matter con- 
tent. In our analysis, while we will allow for both adia- 
batic and isocurvature fluctuations for general Friedmann 
background, we will assume the dark matter content to 
be contributed only by CDM and baryons, ignoring any 
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FIG. 1. Transfer function for the adiabatic (solid curve) 
and isocurvature (dashed curve) modes in the standard cold 
dark matter model (Q m = 1, h = 0.5 and f2 B = 0.05). The 
isocurvature mode is more damped on small scales. 



contribution that may come from neutrinos in the form 
of hot dark matter. 

As for the adiabatic and isocurvature transfer func- 
tions, we take the CDM expressions 

= ln(l + 2.34 g ) 
V ' 2.Mq 

1 + 3.89g + (16. lq) 2 + (5A6qf + (6.71g) 4 
Ti S o(k) = 

(40<?) 2 
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as provided by Bardeen et al. where q = k/T 

( /i _1 Mpc)~ 1 . Here T is the so called "shape" parameter 
and for CDM model is defined as T = hVL m exp(— f2 B — 
VZhQ-s/rim), so as to account for the small-scale damp- 
ing due to the presence of a non-negligible baryon frac- 
tion |25|. Note that we have redefined Eq. jl5| ) with 
respect to the expression given in Ref. so as to have 
the same small-fc asymptotic behavior, T(k) — > 1, for 
both components. Among other things, this means that 
nj SO = 1 corresponds to a scale invariant spectrum, in- 
stead of the usual ni so = —3 appearing in the literature. 
We verified that Eqs. (|l4|) and ( |l5| ) reproduce quite well 
the outputs of the CMBFAST code for the fc-range of 
interest for our analysis. 

As an example, we show in Figure [j] the shape of these 
two transfer functions for a particular choice of the cos- 
mological parameters. It is apparent that the isocurva- 
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FIG. 2. Power spectra of pure isocurvature (dotted curve) 
and adiabatic (dashed curve) modes with spectral indices 
?Jad = 0.962 and rn so = 1.56. The middle (solid curve) power 
spectrum is a mixed one, corresponding to a = 8 x 10~ s . The 
data points are from the APM survey. 



ture case has a larger damping at large k values, due to 
the fact that they start out as zero-curvature or isother- 
mal fluctuations and takes longer than the adiabatic ones 
to build up, after the matter-radiation equality epoch. 

Since adiabatic and isocurvature are independent 
modes of fluctuations, when both modes are present at 
the same time the total power spectra of matter and radi- 
ation can be computed as a linear combination of the pure 
isocurvature and adiabatic ones. In the LM model con- 
sidered in this work, both the adiabatic and isocurvature 
modes have power law initial power spectra, but with dif- 
ferent spectral indices. Therefore, the overall spectrum 
can be casted in the form 

Pmix(fc) = AUx [(kT Q ) n ^T 2 d (k) + a (kT ) n -T 2 o (k)] , 

(16) 

where tq is the conformal time at present and a is a di- 
mensionless parameter that indicates the relative contri- 
bution of isocurvature and adiabatic perturbations (pure 
adiabatic and isocurvature power spectra correspond to 
a = and a — oo respectively) . The comparison of our a 
with the mixing coefficients introduced by other authors 
(e.g. 1 15 l^Jl^]) may not be straightforward because, in 
contrast with their approach, we have also allowed for 
"ad 7^ "iso 7^ 1- In the case n a( j = 7ij SO = 1, our a com- 
pares with the other choices as follows: a = (1 — as) /as 
M|, a = a K KSY @, ol = a E i</(l ~ a E K) In 
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FIG. 3. Mixed CMB spectra with fixed n iso value. The 
solid lines correspond to the pure adiabatic and isocurvature 
modes, normalized to COBE. The dotted-dashed lines cor- 
responds to a = 3, 6, 9, 12 x 10 -3 from top to bottom. The 
points shown here are the binned spectrum from Ref. | p7[ . 
However, we have used the whole set of experiments from 
Table | in the CMB analysis. 

eq.(|l^), iV m i x is the normalization factor, that we de- 
rived by normalizing the CMB spectrum to the COBE 
data |22|,^3]]. Different mixing coefficients a lead to dif- 
ferent normalizations: 



N, 
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(17) 



In the above expression, / = N a d/N iso where iV a( j 
and iVi so are the normalizations of the pure adiabatic 
and isocurvature scalar modes, and their ratio conveys 
the information about the different Sachs-Wolfe con- 
tribution of isocurvature and adiabatic modes to CMB 
anisotropics. In the normalization, we also took into ac- 
count the tensor contribution to the anisotropics, accord- 
ing to the parameter specified below. A typical hybrid 
spectrum is plotted in fig. ||, together with the APM 
data points |2qj . Given the big tilt of the isocurvature 
power spectrum, the mixed one has a slope similar to 
the adiabatic one at small k. The effect of the differ- 
ent normalization is evident even for low values of the a 
coefficient. 

As for the CMB anisotropics, following the standard 
notation, we describe their dependence on the direction 
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FIG. 4. Mixed Ci spectra with different tilts of the isocur- 
vature component with fixed ratio a = 0.012. The solid line 
is the pure adiabatic mode (n a d = 0.962), the dotted line cor- 
responds to mixed isocurvature (a = 0.012) with nj so = 1.36 
and the dashed line to rii so = 1.61. 



Therefore, the radiation power spectrum is defined as 

Q = (\ai m \ 2 ) , (19) 

where the brackets denote the ensemble average over dif- 
ferent realizations. When both scalar and tensor modes 
are present, the Ci can be decomposed as C; = Cf + Cj . 
In Linde-Mukhanov model the tensor contribution is in- 
cluded in the adiabatic mode. 

Similarly to what happens to the matter power spec- 
trum, the radiation power spectrum for the mixed case 
can be found as a linear combination of the two indepen- 
dent adiabatic and isocurvature power spectra Cf d and 

c'r ■■ 



1 + fa 



(20) 



In the expression above Cf and Cf° are both normalized 
to COBE, according to ref. f23|| . In computing the power 
spectra, we modified the CMBFAST code Jl|@ to our 
purposes. In figures |^ and |] we show some examples of 
CMB spectra, together with the estimated binned spec- 
trum from Bond et al. , which provides a visual indi- 
cation of the experimental status. However, we have used 
the whole set of experiments from Table | in the CMB 
analysis. The effect of adding an isocurvature compo- 
nent is to add a lot of power on small multipoles through 
the SW effect. The anti-tilt of the spectrum is in gen- 
eral not enough to compensate this effect, and the first 
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FIG. 5. Typical values of v as a function of M/m, with 
dio in the range 1-10, from top to bottom. 

acoustic peak is consequently lower than in the corre- 
sponding pure adiabatic case. Fig. ^ shows examples of 
mixed spectra with different a values for fixed n iso and 
n a d- It shows that even a small a can cause a signifi- 
cant damping of the acoustic peaks. Moreover note that 
a high n; so component also leads to suppressed acoustic 
peaks, for a fixed a value (see Fig. ||). This is because 
the bigger is nj SO the smaller is the normalization coef- 
ficient iVi so and the bigger is /. Therefore for fixed a, 
the isocurvature Ci spectrum takes more weight in the 
mixture. 



A. Choice of model parameters 

We will consider here a set of values for the parameters 
of the model. Let us start with the scalar spectrum. The 
adiabatic tilt is given by Eq. (|i"o|), 

f»ad = 1 - 2/^hor = 0.9692 , (21) 

for iVhor = 65, while the spectral index for isocurvature 
fluctuations is 

n iso = 0.9722 + 4v/3 . (22) 

Figure || shows some typical values of the parameter v as 
a function of M/m, for different choices of dio- 
As for the tensor contribution, we considered: 

n T = -l/N hoI = -0.0154, (23) 

and a tensor to scalar ratio given by: 

C T 

R=-^~-7n T = 0.1077, (24) 

^2 

which is not negligible. In Fig. ^ the tensor contribution 
is very small and is already included in the solid line for 
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v 

FIG. 6. Typical values of a as a function of the tilt pa- 
rameter v in the isocurvature spectrum. The different curves 
refer to different dio ranging from 1 to 20, from bottom to 
top. 

the adiabatic component, in order to make emphasis on 
its difference with respect to the isocurvature component. 

With the above spectral properties, we determined the 
relative amplitude a in ( |l6| ) to be related to the input 
parameters dio and v, 

a = d w x 10- 1291 - 4 343 " . (25) 

For any value of dio between 1 and 10, the tilt parameter 
v ranges between 0.29 and 0.48. As a consequence, n iso 
is found to be between 1.36 and 1.6lj^| Note that while v 
fixes the value of the isocurvature spectral index, several 
a values are still possible, depending on the dio values. 
In fig. U we plot the value of a as a function of v, for dio 
in the range 1 — 20. In any case, the values of a found 
are small, and for dio < 32, a never exceeds 0.09. In 
the comparison with the data, we considered a < 0.08 
and 1.36 < nj SO < 1.61, and treated them as independent 
parameters, although they are not really independent [see 
Eq. @]. 

Note that we have represented in Fig. g the relative 
contribution a of the isocurvature component to the to- 
tal power spectrum at COBE scales, i.e. large scales. Due 
to the strong positive tilt of the isocurvature spectrum 
(rtiso > 1), this relative contribution increases towards 
smaller scales. Since one of the observational constraints 
that we will consider in the following is represented by the 
number density of local galaxy clusters, it is interesting 
to estimate the isocurvature contribution at the charac- 
teristic cluster scales, ~ 10/i _1 Mpc. To this purpose we 



3 Note that this is equivalent to —2.64 < 7ii so < —2.39 in 
the usual notation, where 7ii so = —3 is the scale invariant 
isocurvature perturbation. 
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FIG. 7. The ratio a c i us /a between the relative isocur- 
vature contribution a at the scale of clusters and at COBE 
scales, as a function of the tilt difference, rii so — n a( j. 



introduce the quantity 

a c ius/a = (A: c i us /fc hor )™ iBO ~ 



(300)"™- 



(26) 



which provides the isocurvature contribution at the scale 
of galaxy clusters. We plot the cluster- to large-scale 
ratio of the isocurvature fraction in Fig. [f| It is apparent 
that the enhancement of the isocurvature contribution 
can be rather large, such that its effect is not negligible 
at the cluster scale, even for a rather small large-scale 
contribution. 

We will discuss in Section [v] the implications of our 
results on the Linde-Mukhanov model. As for the cos- 
mological parameters f2 m , f?A and Hq, we always consid- 
ered flat models (0 TO + = 1), with 0.2 < fi TO < 1, and 
H = 100 h = 50,65,80 km s" 1 Mpc" 1 . As for the den- 
sity parameter contributed by baryons, we take the value 
f2 B = 0.019 h~ 2 , which follows from the low deuterium 
abundance, as determined by Buries & Tytler Pq|. 



IV. OBSERVATIONAL CONSTRAINTS 

In this section we test our model against the available 
microwave background and large scale structure data, 
and we make predictions on the precision with which 
future satellite experiments will determine the relevant 
parameters. To this purpose, we use the Fisher matrix 
technique to determine how well MAP and Planck satel- 
lites will constrain the isocurvature contribution to the 
total power spectrum and its spectral index. 



A. The microwave background data 

In the comparison of the CMB data with the model 
predictions, we performed a x 2 analysis, first applied to 



CMB data by Lineweaver et al. J2j|. More precisely, we 
computed the x' ''■> as a function of the cosmological pa- 
rameters A, on the band-power estimates of the CMB 
data, ST ohs {n), and the model predictions, ST mod {X,n), 



jiven iVexp observed data points with their errors a 



X 2 (A) = £ 



ST ohs (n) - ST mod (X,n) 
cr(n) 



(27) 



Evaluating the expression above, we used the 41 data 
points pp| ] listed in table [|, and the corresponding win- 
dow functions. We chose not to introduce the recent 



points from MAT |3lJ] and Python V 32 experiments 
because the former does not yet include calibration er- 
rors and the results of the latter are still under discussion. 
The value of x 2 obtained is a function of the model pa- 
rameters A. 

For each choice of the parameters VL m and Hq, which 
describe the Friedmann background, we compute the x 2 
between model and data in the (ni so , a) parameter space. 
In Table [n] we report the values of the parameters of the 
best fit, the corresponding Xmin value, and the probabil- 
ity of getting that x 2 value with present data if the model 
considered is the real one. 

For h — 0.5 a small contribution of the isocurvature 
component is desirable, especially for low f2 m universes. 
On the other hand, for h — 0.8 the addition of an isocur- 
vature contribution to the adiabatic mode doesn't pro- 
vide a better fit to the data. In general, however, the al- 
lowed isocurvature fraction tends to be small (a < 0.01), 
while the isocurvature spectral index is not significantly 
constrained (see Fig. ||). In any case, the best fit to the 
data is provided by the lowest 7ii so considered. 



B. Combining CMB and LSS constraints 

In order to further constrain the parameter space of 
allowed models, we will consider in this Section the con- 
straints coming from large scale structure observations. 
In particular, we will constrain the shape of the power 
spectrum, by comparing to results from the analysis of 
galaxy clustering, and its amplitude by resorting to con- 
straints from the local abundance of rich galaxy clusters. 

As for the shape of the galaxy power spectrum, dif- 
ferent determinations have been realized in the last few 
years, both for projected |^6| and redshift |33| samples. 
Such analyses converge to indicate that the observed 
galaxy power spectrum is well reproduced by an adia- 
batic CDM-like P(k), in a flat universe, with shape pa- 
rameter 0.2 < r < 0.3, for a scale invariant primordial 
spectrum Q. 

According to the results of the previous section on the 
CMB constraints, the relative contribution of the isocur- 
vature component of the fluctuations is always rather 
small. As a result, the shape of the purely adiabatic 
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spectrum is never significantly changed by the isocurva- 
ture component. For this reason, in order to implement 
the constraint from the shape of the galaxy power spec- 
trum, we will simply require in the following that the 
adiabatic component of our mixed fluctuation spectrum 
has a shape parameter T lying in the 0.2-0.3 range. 

As for the amplitude of the power spectrum, a powerful 
constraint is represented by the number density of nearby 
galaxy clusters. Since rich galaxy clusters involve a typi- 
cal mass of the order of lO 15 /).^ 1 M©, their number density 
is connected to the amplitude of density fluctuations at 
scales of about 10 /i _1 Mpc. Analytical approaches, based 
on the method originally devised by Press & Schechter 
p5| , show that the cluster abundance actually constrains 
the quantity as — <7gf2£j, where as is the r.m.s. fluc- 
tuation amplitude within a top-hat sphere of 8/i _1 Mpc 
and (3 ~ 0.4-0.5, almost independent of the shape of the 
power spectrum and of the presence of a cosmological 
constant term [^6| . 

Different analyses, based on the distributions of Ar- 
ray cluster luminosities, X-ray temperatures and veloc- 
ity dispersions of member galaxies, converge to values 
of a s in the range 0.5-0.6 |37|. For definiteness, in the 
following we will use for our analysis the expression 

a s = (0.55 ± 0.05) ft-o.43+o.09n m (28) 

where the errors are intended to formally correspond to 
a 90% confidence level, while the £l m dependence is that 
provided by Girardi et al. for a flat Universe, with 

We show in Figure |^ constraints on the nj so — a param- 
eter space by combining results from LSS and CMB data. 
Each panel corresponds to a choice for the (Q m ,h) pa- 
rameters, among those reported in Table which agrees 
with the measured shape of the galaxy power spectrum. 

As for the constraints from the abundance of local clus- 
ters, the 90% confidence level region [see Eq. (|28|)1 is 
shown by the shaded area. We note that, as nj so in- 
creases, the contribution of the isocurvature fluctuations 
at the cluster scale also increases, thus requiring a smaller 
a value to keep the power spectrum amplitude at the level 
required by the cluster number density. 

In order to establish confidence levels for model exclu- 
sion from the analysis of CMB anisotropics we consider 



the quantity A\ 



Xmin ' 



where Xmin * s the min- 



imum value as reported in Table [n| We assume the x 2 
statistics for A^of = 38 to be normally distributed with 
mean Adof and r.m.s. scatter given by V2Ad f ■ Accord- 
ingly, all the models displayed in Figure || correspond to 
an acceptable value of Xmin' The solid curve indicates 
the 90% c.l. upper limit on a, which corresponds to 
A\ 2 = 4.61 for two significant fitting parameters. 

As a general result, it is interesting to note that large 
scale structure constraints significantly contribute to fur- 
ther restrict the range of the allowed parameter space. 
For instance, the two models with (fl m ,h) = (0.5,0.65) 
and (0.4, 0.80) are now ruled out, since they can not sat- 
isfy at the same time both the CMB and local cluster 



abundance constraints, while the model with (£l m ,h) = 
(0.5,0.50) is constrained by the cluster abundance to 
have a < 0.0015. 

It is worth reminding here that the constraint of Eq. 
( p8| ) from the local cluster abundance has been derived in 
the literature under the assumption of Gaussian statis- 
tics for the density fluctuations. On the other hand, the 
density fluctuations predicted by our model are given by 
a scale dependent superposition of a Gaussian adiabatic 
components and of a non-Gaussian isocurvature compo- 
nent, whose probability density function (PDF) corre- 
sponds to a x 2 model with one degree of freedom. Al- 
though the extension of the Press-Schechter formalism 
to non-Gaussian statistics has been pursued by differ- 
ent authors f39j , such attempts concentrated on scale- 
independent PDF models. In our case, we expect that 
the positive skewness of the x 2 -distribution should ease 
the formation of galaxy clusters, for a fixed as, as a conse- 
quence of the broader high density tail in the PDF for the 
isocurvature component. Therefore, the net effect would 
go in the direction of decreasing the required fluctuation 
amplitude at the cluster scale and, thus, to somewhat 
increase the allowed a values. In any case, since the 
isocurvature component is always constrained to be rel- 
atively small even at the cluster scales (a c ius S 0.15; cf. 
Fig. |7|) , we are confident that our assumption of Gaussian 
statistics should be a sensible one. 



C. The future CMB experiments 

In this section we compute the estimates of the er- 
rors with which the future satellite experiments MAP 
and Planck will determine the isocurvature contribution 
to the CMB power spectrum. The aim here is to verify 
whether the future CMB data alone will be able to con- 
strain small isocurvature contributions when other cos- 
mological parameters are constrained at the same time. 

In order to provide such an estimate, we resort to the 
Fisher information matrix approach |40|| . When no po- 
larization is considered, the Fisher information matrix is 
defined as 



„ _ x dC x 



Cov = 



^( c ' + - v, ' i )' 



(21 + l)/ sky 



(29) 
(30) 



Here / s k y is the fraction of sky covered, <7f, is the Gaus- 
sian beamwidth (a\, — # fwhm /V81n2, and # twhm is the 
full width at half maximum), \ is the set of parameter 
whose errors have to be determined, and w contains the 
information on the detector resolution and sensitivity: 



where 



pixel 



a pixcl^P ixel 
1 



(31) 
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FIG. 8. Constraints on the rii so — a plane from cluster abundance and CMB anisotropies. Each panel corresponds to a 
choice for the (f2 m , h) parameters which satisfies the constraints on the shape parameter F. The corresponding F values for 
each model are also reported. The shaded area correspond to the 90% c.l. from the cluster abundance, obtained according to 
eq.(ps|). The solid curves indicate the 90% c.l. upper limit on a from the CMB constraints (see text); oscillations in the shape 
of these curves are due to limitations in the numerical precision. 



In the case that also polarization is considered, the 
expression for the Fisher and covariance matrices be- 
come more complicated, and we refer here to Ref. pd] ] 
for the explicit expression. We note that in the ex- 
pression above a perfect foreground subtraction is as- 
sumed, so that the estimates found should be considered 
somehow ideal. Also notice that the expression reported 
here strictly holds for Gaussian perturbations. Therefore, 
its application to the Linde-Mukhanov model should be 
taken with some caution, especially in the very large mul- 
tipole regime (i.e., small scales), where the non-Gaussian 
isocurvature contribution may be non negligible. 

We estimate the expected errors on a and ni SO with and 
without polarization. In the calculation of the Fisher 
matrix, we considered as free parameters the Hubble 
constant h, the baryon abundance Ct B h 2 , the normaliza- 
tion Cio, the reionization optical depth r, the isocurva- 
ture/adiabatic ratio a, and the isocurvature spectral in- 
dex rii S0 . Eqs. ( |30| ) and ( p9|) are estimated by computing 
Ci for a fiducial model with O m = 0.3 (fixed), h = 0.65, 
B /i 2 = 0.019, Cio = Ci 0cOBE and r = 0.05. 

The precision in the determination of a and nj so is 
reported in table |^ for different fiducial values, and dif- 
ferent experiments. For both experiments, we combine 
different channels in order to give the estimates on the er- 
rors. More precisely, we consider the 3 highest frequency 
channels for MAP (40, 60, 90 GHz), and the 5 central fre- 
quency channels from Planck (70, 100 GHz channels from 
LFI and the 100, 143, 217 GHz channels from HFI). The 



values of the beamwidths and sensitivities for the two 



experiments are reported in table III, according to Refs. 
il- 

As a result, we find that the future satellite experi- 
ments will be able to constrain the spectral index rii SO 
much better than the present data. Clearly, nj SO is bet- 
ter determined if a higher isocurvature component is al- 
lowed, and for a given a a higher fiducial isocurvature 
anti-tilt ensures a smaller uncertainty in the parameter 
determination. Considering polarization is also useful in 
this respect, and typically reduces the errors by almost 
a factor of 2. If the isocurvature contribution is as low 
as indicated by our present analysis, rii SO will be anyway 
determined by Planck with an error of about 0.1. 

As for the constraints on a, table IV shows that if the 
isocurvature contribution is quite large (a ~ 0.015) both 
MAP and Planck will be able to detect it; on the con- 
trary, if a is as low as the present data seem to suggest, 
only Planck (with polarization information included) will 
be able to claim a detection, while MAP will only put an 
upper limit at the 0.002 — 0.003 level. Note that the use 
of polarization data certainly helps in reducing the error, 
although we don't find the big improvement claimed by 
p8[ . This may occur because we keep the scalar to ten- 
sor ratio R fixed in this model, and therefore we are not 
affected by the degeneracy of this parameter with a. Po- 
larization is useful in breaking the degeneracy between 
the reionization parameter r and the isocurvature con- 
tribution a and helps in reducing the errors especially in 







the case of the Planck experiment. 



V. DISCUSSION 



knowledge SISSA for hosting all of them during different 
phases of this work. The authors thank Andrei Linde, 
Andrew Liddle, and Douglas Scott for generous discus- 
sions. 



In this paper we have investigated the consequences 
of a CDM cosmology with mixed isocurvature and adia- 
batic initial conditions as prescribed by a generic Linde 
Mukhanov inflationary model. 

We showed how the total spectra are modified by the 
isocurvature contribution, as both the primordial spec- 
tral index for the isocurvature power spectrum, nj so , and 
the isocurvature fraction, a, are varied within their al- 
lowed ranges. In order to constrain the parameters of 
the model we resorted to available data on the CMB 
anisotropy, as well as on the large-scale structure con- 
straints from the shape of the galaxy power spectrum 
and the number density of nearby galaxy clusters. Ob- 
servational constraints from CMB and LSS data have 
been shown to provide complementary informations. As 
a result, we found that the allowed isocurvature contri- 
bution at COBE scales is always very small, a < 0.006. 
Therefore, our results generalizes to the nj so > 1 case and 
strengthen the conclusions reached in Ref. [jl6| based on 
LSS constraints alone. 

We note that, even allowing for the strong positive 
tilt of the isocurvature component, the permitted isocur- 
vature contribution at the cluster scales is always small, 
with a c i us < 0.15. However, this contribution is expected 
to increase to ~ 50% level when we consider smaller 
scales, < 1 /i _1 Mpc, which are relevant for galaxy for- 
mation. Therefore, the resulting x 2 (positive skewness) 
non-Gaussian statistics contributed by the isocurvature 
fluctuations can play a significant role to ease the galaxy 
formation at high redshift. The phcnomcnological impli- 
cations on galaxy formation of a small isocurvature con- 
tribution at the COBE scales remains to be investigated 
in detail. 

As for the determination of the isocurvature spectral 
index nj so , the best fit to the CMB data is always pro- 
vided by the lowest nj so considered; although the limits 
from LSS and CMB data in Fig. || show a mild depen- 
dence on n- lso . In this respect, future CMB experiments 
can help to set more accurate limits on rii so , with an op- 
timistic estimate of 3% and 1% (high ni so and a) and a 
more realistic one of 20% and 8% for MAP and Planck 
respectively. 

Finally, we showed that only the Planck experiment 
will be sensitive enough to detect a possible non vanishing 
value for the a parameter within the limits already set 
by present CMB and LSS data. 
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OO.U_ 22 g 


53 


Femenia et al (1998) 


IAB 


04 5+41-8 


125 


Piccirillo & Calisse (93) 


MSAM 


62.0+2l' 7 


143 


Cheng et al (1996) 


MS AM 


60.4i 2 2 ^ 


249 


Cheng et al (1996) 


MSAM 


5o.oli; : g 


160 


Cheng et al (1997) 


MSAM 


65.0111;° 


270 


Cheng et al (1997) 


QMAP 


47.01°' 


80 


deOliveira-Costa et al (98) 


QMAP 


59.01°; 


126 


deOliveira-Costa et al (98) 


QMAP 


52.0l^ : 


111 


deOliveira-Costa et al (98) 


OVRO 


56.01°;* 


589 


Leitch et al (1998) 


TABLE I. 


Data points used 


in the x analysis. First col- 



umn is the experiment, second column is the experimental 
value with its error, third column is the effective multipole 
number, and fourth column is the reference paper. 







a 


Xmin 


P(X < Xmin) 


0.2 


1.359 


0.006 


21.6 


0.015 


0.3 


1.359 


0.003 


21.9 


0.017 


0.4 


1.359 


0.001 


22.7 


0.023 


0.5 


1.359 





23.5 


0.032 


0.6 


1.359 





24.8 


0.048 


0.7 


1.359 





26.4 


0.078 


0.8 


1.359 





28.3 


0.126 


0.9 


1.359 





30.5 


0.200 


1 


1.359 





33.5 


0.323 


0.2 


1.359 


0.001 


24.4 


0.043 


0.3 


1.359 





26.4 


0.08 


0.4 


1.359 





29.2 


0.15 


0.5 


1.359 





32.7 


0.29 


0.6 


1.359 





37.4 


0.50 


0.7 


1.359 





43.1 


0.74 


0.8 


1.359 





49.2 


0.89 


0.9 


1.359 





55.5 


0.97 


1 


1.359 





61.8 


0.99 


0.2 


1.359 





30.2 


0.19 


0.3 


1.359 





36.0 


0.44 


0.4 


1.359 





44.1 


0.77 


0.5 


1.359 





52.8 


0.94 


0.6 


1.359 





61.5 


0.99 


0.7 


1.359 





69.8 


0.99 


0.8 


1.359 





77.8 


0.99 


0.9 


1.359 





85.4 


0.99 


1 


1.359 





92.6 


0.99 



TABLE II. Results of the \ 2 analysis for different models 
(h = 0.50,0.65,0.80 from top to bottom). Column 1: the 
matter density parameter. Column 2 and 3: the best fit values 
for rii so and a respectively. Column 4: the corresponding 
Xmin value. Column 5: the probability of getting a smaller 
X 2 , assuming that this is the real model. We considered here 
41 experimental points and 3 free parameters (a, n,i so , and 



the normalization), which corresponds to iVdof = 38 degrees 
of freedom. 



Experiment frequency (GHz) # twhm cr p (fjK) Op° l (/iK) 

MAP 40 0.47° 35 49.3 

MAP 60 0.35° 35 49.3 

MAP 90 0.21° 35 49.3 

Planck-LFI 70 14' 9.8 13.8 

Planck-LFI 100 10' 11.72 16.5 

Planck-HFI 100 10.7' 4.63 

Planck-HFI 143 8.0' 5.45 10.2 

Planck-HFI 217 5.5' 11.7 26.2 



TABLE III. Beamwidths and sensitivities of the satellites 
experiments MAP and Planck. 
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fiducial 


fiducial 


MAP 


MAP 


Planck 


Planck 


a 




5a 


<5?liso 


5a 


5ni BO 


0.002 


1.439 


0.0035 


0.315 


0.0024 


0.209 


0.002 


1.572 


0.0022 


0.195 


0.0016 


0.139 


0.008 


1.439 


0.0044 


0.094 


0.0031 


0.064 


0.008 


1.572 


0.0030 


0.062 


0.0023 


0.047 


0.015 


1.439 


0.0054 


0.059 


0.0040 


0.042 


0.015 


1.572 


0.0039 


0.041 


0.0031 


0.031 


0.002 


1.439 


0.0025 


0.297 


0.0013 


0.126 


0.002 


1.572 


0.0015 


0.172 


0.0008 


0.073 


0.008 


1.439 


0.0031 


0.091 


0.0017 


0.039 


0.008 


1.572 


0.0020 


0.056 


0.0010 


0.024 


0.015 


1.439 


0.0038 


0.059 


0.0020 


0.025 


0.015 


1.572 


0.0026 


0.038 


0.0014 


0.016 



TABLE IV. Estimation of the isocurvature mode parame- 
ters with the future satellite experiments MAP and Planck. 
Top: polarization not considered. Bottom: polarization con- 
sidered. The first two columns indicate the fiducial value 
considered, and the others indicate the estimated errors in 
the parameters. 
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